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ABSTRACT 

A unified model is proposed for the radio and X-ray outburst of nova CI 
Cam 1998 which suggests the shock interaction of nova shell with the circum- 
stellar gas. The spherical model is able to describe kinematics of the radio shell 
together with the evolution of the radio and X-ray fluxes. However, the X-ray 
spectrum in this model is harder than the observed one. Better agreement with 
observations demonstrates the model in which the spherical shell interacts with 
the nonspherical circumstellar medium. The latter is made up of the broad bipo- 
lar jets with the openning angle of 120° and the dense equatorial wind. In the 
optimal model the kinetic energy of the nova shell is ~ 8 x 10 43 erg, while the 
shell mass lies in the range of (1 — 5) x 1CT 7 M & . 



Introduction 



Nova CI Cam was discovered 1998 March 31.64 by the ASM det ector of RXTE as 



an X-ray outburst wit h the flux of 139 mCrab in the range 2-12 keV (jSmith et al.l 11998 



Revnivtsev et al.lll999l ). At the earlier epoch March 31.36 this source was absent with the 
upper limit of 40 mCrab. The extrapolation of the flux evolution backwards results in the 
outburst moment March 31.5, i.e., 50903.5 MJD; this is considered below as zero point of 
the nova shell expansion. The X-ray s ource was identified with the binary of symbiotic 



type CI Cam (IWegner &: Starrfieldlll998l). in w hich the luminosity is dominated by B-giant 



classified as B4 III-V[e] (IBarsukova et al. 
eccent ricity e 



20061 ) . The binary has the orbital period of 19.4 d, 
0.62 and maj or semiaxis a = 4.8 x 10 12 / sin i cm, where i is the inclination 



angle (IBarsukova et al.l 120061 ). The optical outburst was observed starting from April 3 at 
the decline stage with the exponential time of 3.4 d (IGarcia et al.lll998t IClark et al.l 120001 ) . 



Nova CI Cam 1998 turned out to be a bright s ource of the radio emission that was ob- 
served on several frequences starting from April 1.9 ( lHjelmingill998l ; iMioduszewski fc Rupen 
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2004 ). The source mapping in radio band during several consequitive epochs staring from 
the beginning sho ws the source expansion with the velocity of ~ 12000 km s" 1 assuming the 
distance of 5 kpc (IMioduszewski &: Rupenll2004l ). The radio emission is interpreted as a syn- 
chrotron emission of the shock wave propagating in the circumstellar medium (CSM). The 
radio evolution shows the strong early absorption at 2 GHz which c an be caused by the syn- 
chrot ron self-absorption (SSA) or/and //-absorption in the CS gas (IMioduszewski fc Rupen 
2004 ). At the late stage t ~ 80 d the radio shell shows an oval shape which is inter- 
preted as a bipolar structure which can be interpreted as a prolate spheroid, bipolar jets, 
or bipolar spherical c louds (IMioduszewski fc Rupen! 12004 ). Interestingly, in this regard, 
(IFilippova et al.l 120081 ) note that the spherical shock wave in the wind predicts significantly 
slower decay of the the X-ray emission compared with observations and that the better 
agreement shows the model of the shock wave in the equatorial disk. 

The nature of the unseen stellar component of CI Cam and the origin of the outburst are 
unclear. There is a conjecture that n ova CI Cam is the result of the outburs t in the binary 
containing black hole or neutron star (IFrontera et al.lll998l ; iBelloni et al.lll999l ). On the other 
hand there are serious arguments that unseen component is the white dwarf (WD) and that 
the nova CI Cam 1998 is the res ult of the thermonuclear explosion on the surface of the 



white dwarf as in classical novae (llshida et al.l 12004 ; iBarsukova et al.l 120061 ; iFilippova et al. 
2008T) . 



Keeping in mind unusual nature of nova CI Cam it is desirable to advance in the 
interpretation of data on the radio and X-ray outbursts. Particularly it would be of interest 
to construct a model of the interaction of the nova shell with the CS gas that might account 
for the radio and X-ray data simultaneously. Such unified model possibly might provide 
us with interesting estimates of the crucial parameters, e.g., energy and mass of the nova 
shell. The radio and X-ray data indicate the possibility of the nonspherical distribution of 
the CSM. It would be sensible, therefore, to explore both spherical and nonspherical version 
of the model. This paper presents an attempt of the construction of the unified model in 
the framework of spherical and nonspherical approaches. 

Model parameters essentially depend on the adopted distance (d) which is poorely 
known. Most p opular are two options: d = 2 kpc dClark et al.ll2000l ; IFilippova et al.l 120081 ) 
and d = 5 kpc (IMioduszewski fc Rupenl |2004 ; llshida et al. 120041 I will present results for 
both distance choices. 
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2. Dynamical model 



The basic model suggests the shock interaction of a spherical nova shell with a spherical 
CS medium. This model with some modification can be also used in the nonspherical case. 
Because the mechanism of the initial energy release is not completely clear, there is some 
uncertainty in the initial density and velocity distribution along the radius in the nova shell. 
Two mechanisms of the shell ejection are conceivable: an explosion which suggests the energy 
release on time scale less or comparable to hydrodynamic time scale and a relatively slow 
matter acceleration by the pressure gradient (optically thick wind). The latter is a major 
mechanism for the mass loss in classic novae (jPrialniki Il986l ; iKato fc Hachisul 119941 ) . In case 
of CI Cam this mechanism however cannot explain the high expansion velocity of the radi o 
source which exceeds 5000 km s _1 for the distance > 2 kpc (IMioduszewski fc Rupenl 12004 ). 
I therefore assume that the nova shell was ejected explosively. In this case one expects that 
the shell expands homologously soon after the passage of the shock wave, i.e., in the time of 
~ 10 2 s. However, this does not specify the density distribution. Moreover, we do not know 
the matter distribution in the close vicinity of the white dwarf that could modify the initial 
density distribution of the nova shell. 

Bearing in mind above uncertainties it is reasonable to use a most simple dynamic model 
in order to minimize the number of free parameters; the Sedov point explosion model is most 
appropriate. However one has to take into account also the initial deceleration phase. Indeed, 
even for the low mass ejecta ~ 3 x 10~ 7 M & expanding in the dense CSM n ~ 10 10 cm -3 
the deceleration length scale is large i?~2x 10 13 cm with the corresponding time scale of 
~ R/v ~ 1 d. The deceleration phase is taken into account assuming the nova shell brakes 
as a whole with the energy conservation. Equations of the energy and mass conservation is 
solved by Runge-Kutta forth order method. 

Note, the formulated deceleration problem for the shell of mass Mq and the energy 
E = (1/2)MqVq expanding in the stationary wind with the density p = w/4irr 2 can be 
solved analitically 

t \ 2/3 



'1 



where r s = Mq/w is the radius at which the mass of the swept up shell becomes equal to 
the initial shell mass, and t s = (2/3)r s /v . At times t 3> t s this solution approaches Sedov 
solution r oc t 2//3 . For Mq = 3 x 10~ 7 M , and wind density parameter w = 10 13 g cm -1 
the characteristic radius is r s = 6 x 10 13 cm. For the initial velocity vq = 6000 km s" 1 
the expansion enters the Sedov stage after about one day. The equation ([1]) is used to test 
numerical solution. 



-4- 



Following iFilippova et al.l (120081 ) I adopt the CSM distribution presented by a homo- 
geneous core and an outer stationary wind. Specifically, p = Po = const in the inner zone 
r<r = 2xl0 13 cm and p = Po(r/r )~ 2 for r > r . The homogeneous core approximately 
takes into account the orbital motion of the less massive component with the velocity of ~ 230 



km s 1 which affects the slow wind flow with the velocity of u = 32 km s 1 (jRobinson et al. 



20021 ) in the orbital zone r ~ 10 13 cm. In fact, computations show that the shell dynamics 



only weakly depends on the variations of the density distribution in the region r < r . 

A large volume of the parameter space was explored; below only two versions of the 
spherically-symmetric model are presented. They correpond to the minimal energy for which 
we are able to describe the expansion kinematics of the radio shell and to reproduce maximal 
X-ray flux for the adopted distance. The shell mass is poorely constrained because the 
expansion dynamics and X-ray ray emission are sensitive primarily to the energy and the CS 
density, not to the shell mass. In the model A for d = 2 kpc and the energy 1.1 x 10 44 erg the 
shell mass lies in the range (1 — 5) x 10~ 7 M Q with the accepted optimal value of 2 x 10 -7 M 
defined by the radius of the radio source at the age t ~ 2 d. In the model B for d = 5 kpc 
with E — 1.6 X 10 45 erg I adopt the same mass 2 x 10~ 7 M . Models parameters are listed in 
Table 1 which, starting from the second column, shows the distance, shell mass and energy, 
CS density in the central zone, r < r , parameter of the density of the relativistic energy e r , 
spectral index of the relativistic electrons, and the filling factor of the radio-emitting clouds. 
The sense of some parameters is desribed in the appropriate place below. 

The evolution of the shock wave radius for the models A and B (Fig. 1 and Table 
1) satisfactorily reproduces the observed expansion of the radio source at the time scale of 
~ 10 2 d. The initial velocity in the model A is 7400 km s _1 , somewhat larger than maximal 
velocities of novae shells but st ill within escape ve locities ~ 1.1 x 10 4 km s~ x for white 
dwarfs with masses ~ 1.3 M & (jAlthaus et al.l 120051 ). In the model B with d — 5 kpc the 
initial velocity 2.8 x 10 4 km s _1 is significantly larger than the velocity expected for any 
known type of cataclismic system with the white dwarf. 



3. X-ray emission of the shock wave 

In our model the reverse shock is absent. Fortunately, the role of the reverse shock in 
the X-ray emission is negligible: it is important only at the early decelerat ion stage t ~ 1 



d. M oreover even at this stage it contributes only in the soft band < 1 keV ( iFilippova et al. 



20091 ) because the speed of the reverse shock is relatively low. The X-ray emission of the 



forward shock is calculated assuming a homogeneous density of the shocked gas. The shock 
speed is assumed to be equal to the shell velocity found from the dynamics computation. 
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The post-shock density is determined assuming strong shock limit. The contribution of the 
relativistic pressure (relativistic particles plus magnetic field) with a fraction 77 o f the total 
press ure is taken into account. The compression factor in this case is 4 + 3rj ( IChevalier 
19831 ). The inclusion of the relativistic pressure requires the itterative procedure: at the first 
step we compute dynamics and X-ray emission; afterwards we calculate radio emission which 
provides us with the estimate of rj and then again recalculate X-ray emission. The electron 
temperatures is calculated from the equation of the heat exchange between ions and electrons 
dT e /dt = Ti/t eq} where t eq = CT^ 2 / rtj is equilibration time, is the ion concentration, C 



is the factor specified in (jSpitzerl Il962l ). The initial temperatures in the shock T e and Tj 
are taken according to Rankine-Hugoniot jump conditions w hile the age is adopted as a n 
integration interval. The cooling function is taken according to ISutherland fc Dopital ( 119931 ). 



The X-ray flux density is presented by the spectrum F(E) oc E exp (-E/T) (where E 
and T are in keV). Calculations show that always T e < Tj. 

The model unabsorbed flux in the range of 3 — 20 keV is shown in F ig. 2 together 
with the unabsorbed flux derived from observations (IFifippova et al.l 120081 ) . Insets show 
the evolution of the co mputed T P compared to the electron temperature recovered from 
the observed spectrum ( Filippova et al.ll2008l ). The flux in the model A is consistent with 
observations only around maximum during the first day; afterwards the flux decays more 
slowly than the obser ved one. The similar behavior shows the model flux in the paper 
( IFilippova et al.l 120081 ). The temperature significantly exceeds the observed one (Fig. 2) 
which reflects higher shock speed than is needed to account for the X-ray spectrum of CI 
Cam 1998. In the model B we face the same problem of the slow flux evolution as in model 
A, while the temperature shows even larger deviation from the observed one compared to 
model A. This difference is related with the substantially larger expansion velocity in the 
model B. Given disparity makes large distance (d > 2 kpc) highly unrealistic. 



4. Modelling radio emission 



The spectrum and high brightness temperature of the radio emission of nova CI Cam 
1998 show that the radiation has a synchrotron origin, while the structure and kinematics 
of the radio so urce idicates that the emissio n originates from the shock wave propagating 
in the CS gas (IMioduszewski fc Rupenll2004j ). The theory of the electron acceleration and 
generation of magnetic field in the shock wave cannot confidently predict parameters that 
define the radio emission. I therefore use a parametric description. 



Below we assume that the magnetic field and relativistic electrons are homogeneously 
distributed on average behind the shock front within the layer of h = 0.2R where R is the 
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shock radius. The equipartition between nagnetic field and cosmic rays with the electron- 
to-proton ratio K ep = 0.01 is adopted. The ratio of the relativistic energy (magnetic field 
plus cosmic rays) to the internal energy in the shock wave (e r ) is a free pararmeter. The 
power spectrum of relativistic electrons in the energy range E > E min = 0.5 MeV is assumed 
dN/dE oc E~ p , where spectral index p is a free parameter. The model takes into account 
the SSA, the //-absorption in the fully ionized wind with the temperature 10 4 K; the Razin 
effect is also taken into account. The latter results in the exponential suppression of the 
spectrum at low frequencies v < u R , where z/ R = (2/3) v^Ivb {y v is the plasma frequency, v B 
is t he gyrofrequ ency) . The Razin effect can be taken into account using the flux prefactor 
(cf. ISimorJll969h 

G = C(z,)exp[-(3/2) 1 / 2 z, R /H, (2) 

where C(v) = 1 for v > z/ R and C(v) = (z^/Vr)^ 2- ^ f° r v < U R- The modelling shows that 
the //-absorption and Razin effect are of little importance compared to the SSA. 

Preceding results, let us show that a homogeneous distribution of the radio-emitting 
material in the shell is inconsistent with the observational data. It is well known that the 
angular size or linear size (if t he distance is available) of the opaque radio source can be 
well determined from the flux (jSlishlll963l ). The size thus obtained weakly depends on the 
magnetic field. Assuming the equipartition one can exclude the magnetic field; moreover the 
size turns o ut to be insensit ive to the ratio of magnetic field and relativistic electrons energy 
density (cf. IChevalierl 1 19981 ) . The dependence of the synchrotron luminosity at 2.25 GHz of 
the expanding source in the luminosity maximum on the peak time t m is shown in Fig. 3 . 
This diagram is similar to that used for the analysis of radio supernovae (jChevalierl 119981 ) . 
The plot demonstrates that for the shell with the homogeneous brightness the observed 
monochromatic luminosiy of CI Cam 1998 in the maximum on day 5 at the distance 2 < d < 5 
kpc corresponds to the expansion velocity 700 — 1600 km s _1 , significantly lower compared 
to the observed value ~ 4000 — 6000 km s" 1 at this epoch (Fig. la). This disparity means 
that the brightness distribution should be essentially inhomogeneous, i.e. the source has a 
cloudy structure. Note that th e inhomogeneity of the radio brightness of CI Cam 1998 is 
apparent from the radio image (IMioduszewski fc Rupenl 120041 ) . 



The description of the inhomogeneous radio source requires additional parameters: fill- 
ing factor (/) and the ratio of the cloud radius to the width of the radio-emitting shell 
(£ = a/h). The intensity of the emergent radiation is then determined by the expressions 



I v = S u [l - exp (-r eS )) , r cS = r g [l - exp (-r c )] 



(3) 



where S = j u /k u is the source function, r g = (3/4)(//£) is the occultation optical depth 
(i.e., the average cloud number along the radius), r c = (4/3)0^ is the average optical depth 
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of the cloud. The surface fraction covered by clouds (covering factor) is 1 — exp(— r g ). The 
expression ([3]) is used for the flux computation, which is sensible approximation. 

The sensitivity of the radio flux behavior to the major parameters is shown in Fig. 4. 
The reference model A is characterized by e r = 0.45, / = 0.025, and £ = a/h = 0.5; three 
other cases show models in which each parameter is twice as smaller. The reduction of e r 
decreases the flux at the optically thin stage and shifts the maximum towards the earlier 
epoch; the reduction of / decreases the flux in all frequences; the reduction of £ shifts the 
maximum towards the earlier epoch since the cloud optical depth decreases. In the latter 
case the flux at the opaque stage is higher because the covering factor gets larger. This plot 
shows how one can recover the optimal model via parameter variations. 

The simulation results for the models A and B are shown in Fig. 5, while the parameter 
values are given in Table 1. In both models the relative cloud radius is assumed to be £ = 0.5 
in order to secure maximal SSA for the moderate parameter e r . Even with this prerequisite 
the required fraction of relativistic component is rather high e r = 0.45 and 0.23 in models 
A and B respectively. As a result, the shock compression factor turnes out to be > 4 



so the fraction of the thermal pressure becomes notably lower than unity (IDrury fc Volk 



19811 ; IChevalierl Il983l ); this fact is taken into account for the calculation of the ion and 
electron temperatures. The filling factor / fa 0.025 in both models while the covering factor 
~ 0.1 (front and rear sides of the shell are included). Models satisfactorily reproduce the 
flux evolution at different frequences, although at late stage the flux decreases more rapidly 
compared to the observations. Maximum at 2.25 GHz is fully determined by the SSA whereas 
the optical depth due to //-absorption is only 0.15 at t = 1 d. The brightness temperature 
of radio-emitting clouds at 2.25 GHz in the flux maximum is ~ 2 x 10 10 K, markedly lower 
than the Compton limit ~ 10 12 K. 

The relativistic energy fraction e r in the preferred model A is rather high (0.45). This 
means that the kinetic energy in this model is close to the minimal permissible value. On 
the other hand, the energy of the model A cannot be significantly larger because otherwise 
we would come to the large X-ray luminosity. The recovered parameters of the model A 
therefore are close to the optimal values for the spherical models. 



5. Nonspherical model 



The bipolar structure of the CI Cam 1998 radio map at late times (IMioduszewski fe Rupen 

20041 ) indicates that CSM might have two-component structure, e.g., rarefied medium in po- 
lar directions combined with the dense equatorial disk. In this picture the radio emission 
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could be associated with the fast shock waves running in the polar directions, while X-rays are 



presum ably emitted by slow shocks propagating in the dense equatorial disk. iFilippova et al. 



( 120081 ) already proposed the equatorial disk for the interpretation of X-ray emission of CI 
Cam 1998. This model is preferred compared to the spherical one because it is able to pro- 
vide the fast decay of the X-ray flux. The ready scenario for the formation of the proposed 
structure of the CSM is absent. However, one could think that the combination of the slow 
wind with the velocity of 32 km s _1 and of the fast disk wind with the velocity of ~ 10 3 km 
s _1 taken together with the gravitation perturbatio n from the unseen component with the 



orbital velocity of 230 km s 1 (IBarsukova et al.ll2006l ) might maintain the suggested structure 
of the CSM. 

To simulate the radio and X-ray emission we consider the CSM composed of the rarefied 
bipolar outflow in cones with the opening angle 29q and dense equatorial disk with the polar 
angle range 9q < 9 < tt — 9$, where 9 is counted form the polar axis z. The density in 
polar cones is constant (po) for r < tq = 2 x 10 13 cm and p = po( r o/ r ) 2 for r > Tq. In 
the equatorial disk the radial density distribution is similar but values are larger. The disk 
height over equatorial plane is assumed to be z oc r for r < r and z = r cos 9 = constant 
for r >> r . This behavior is approximated by the expression 

rr cos O , A s 

z o = • 4 

r + r 

We assume that the spherical nova shell interacts with the asymmetric CSM described above. 
The radio and X-ray emission is calculated separately for the polar cones (model Cp) and for 
the equatorial disk (model Ce) emploing the spherically-symmetric model. The total flux is 
a superposition of both components. The resulting average electron temperature of X-ray- 
emitting plasma is defined as the flux weighted mean of temperatures of both components. 
It should be noted, however, that the contribution of X-rays from the bipolar component 
is small, so the flux and temperature is dominated by the disk shock. Assuming that the 
expansion of the radio source corresponds to the bipolar shocks we must admit that the 
inclination angle obeys to the condition i > 90° — 9q. 

The nonspherical model describes the X-ray flux and the electron temperature (Fig. 6) 
significantly better compared to the spherical model A (Fig. 2). Parameters of the CSM for 
the bipolar cones and the disk are given in Table 1. The X-ray emission is dominated by the 
disk shock, while the radio is related primarily with the bipolar cones. The relativistic energy 
fraction is large (e r = 0.9), so the model energy (8 x 10 43 erg) is close to the minimal value; 
the initial velocity of the nova shell in the model is 6300 km s _1 . Interestingl y, the recovered 



energy is only 1.7 times larger than the energy of the spherical model in IFilippova et al. 



( 120081 ) for the case of the Sedov strong explosion models. 
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Nonspherical model, alleviating the disparity between the model and observed rates of 
the flux decay, does not remov e completely th i s cont radiction. The similar contradiction 



is present in the disk model of iFilippova et all (120081 ) . The fast decay of the observed X- 



ray flux may be related with another effect being absent in the disk shock model, namely, 
the spherization of the post-shock flow due to outflow of shocked gas perpendicular to the 
equatorial plane. However, to take into account this effect requires the multi-dimensional 
hydrodynamics. 

The study of possible uncertainties implies that the nova energy and the density of 
the CSM are determined in our model with the accuracy of about 30%. The mass of the 
nova shell is less certain because results weakly depend on the mass for the fixed value of 
energy. The optimal mass range (1 — 5) x 10 -7 M & is determined from the requirement of the 
description of the early (t ~ 2) expansion. Finaly, it should be stressed that the proposed 
nonspherical model is illustrative; the actual distribution of the CSM may somewhat differ 
from the recovered optimal model. 



6. Discussion 

Our goal was the possibility to account for the radio and X-ray emission from nova CI 
Cam 1998 in the framework of a unified model of the nova shell interaction with the CSM. 
I show that the interaction of the spherical nova shell with the spherical CSM is able to 
describe the the evolution of the radio flux and the expansion kinematics of the radio source 
and only roughly to describe the X-ray emission assuming the distance of 2 kpc. The larger 
distance is less likely because it suggests the larger expansion velocity of the radio source 
and thus more harder X-ray spectrum in odd with observations. For the distance of 2 kpc, 
however, more reasonable description can be found in the model of the interaction of the 
spherical shell with the nonspherical CSM. The latter is composed by the combination of 
the rarefied bibolar outflow and the dense equatorial disk. 

The rarefied gas in bipolar cones can be associat ed with the fast bipola r wind charac- 



teristic of the binary systems with accreation disks. iRobinson et al.l (120021 ) report on the 
detection in the HST spectra taken in March 2000 of the ultraviolet lines of C IV and Si IV 
with P Cygni profiles and the expansion velocities of ~ —1000 km s" 1 in the blue absorption 
wing. With this velocity and wind density from Table 1 the mass loss rate related with the 
bipolar outflow is the n ~ 2.4 x 10 -6 M yr" 1 . For the equatorial wind velocity 32 km s _1 



( IRobinson et al.ll2002[ ) the mass loss rate related with the equatorial wind turnes out to be 
~ 0.5 x 10 ^ Mq yr ^. The total mass loss rate is therefore ~3x 10~ 6 M Q yi- 1 , close to 
~ (1 — 2) x 10 -6 M Q yr _1 obtained for the spherical model of X-ray emission of CI Cam 
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1998 (IFilippova et all 120081 ). 



Given the openning angle of bipolar outflow ~ 120° the inclination angle of the binary 
system should be i > 30° to meet requirement that the shock tangential velocity should 
coinside with the physical expansion velocity. This inequality is consistent with another 
inequality i > 38° found earlier from the requir ement that the unseen component of the 
binary is the white dwarf (IBarsukova et al.ll2006l ). 



Parameters of the optimal model of the nova CI Cam 1998 £~ 8 x 10 43 M ~ (1 — 5) x 
10~ 7 M Q , bipolar structure, and the expansion velocity ~ 6300 km s^ 1 are reminescent of 



the similar properties of recurrent nova RS Oph 2006. Recently iRupen et al.l (120081 ) studied 
bipolar structure of the radio source of RS Oph 2006 and found that the expansion velocity 
of bipolar shocks is ~ 10 4 km s _1 . They estimate the energy of bipolar outflow of RS Oph 
2006 to be ~ 10 44 erg, which is close to the energy of CI Cam 1998; the estimated ejecta 
mass ~ 10~ 7 M Q of RS Oph 2006 is also close to the shell mass of CI Cam 1998. 

Yet there are some apparent differences between these novae. For comparable distances 
(~ 2) kpc and comparable extinction (Ay ~ 2.5 mag) the optical flux in V band at maximum 
of CI Cam 1998 is ~ 40 times lower compared to RS Oph 2006. On the other hand the X-ray 
flux of CI Cam 1998 in 0.5-10 keV band is ~ 20 times larger, while radio flux at 5 GHz is 
~ 7 larger compared to RS Oph 2006. Binary systems are dissimilar as w ell. For CI Cam 



the b inary period is 19 d and the donor is B-giant with the mass > 12 M Q (IBarsukova et al. 



20061 ). whereas for RS Oph the binary perio d is 460 d and the donor is M-giant with the 
mass ~ 0.5 M (IDobrzycka fc Kenyonlll994l ). All these differences leave open the issue of 
the possible simiarity of the outburst physics of CI Cam 1998 and RS Oph 2006. 
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Table 1: Model parameters 

Model d M E p e r p f 

kpc 10" 7 M 10 44 erg 10~ 15 g cm- 3 

A 2 2 LI 2 0.45 2.3 0.025 

B 5 2 16 3.2 0.23 2.3 0.025 

Cp 2 2 0.8 0.6 0.9 2.1 0.04 

Ce 2 2 0.8 4 0.27 2.1 0.04 
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Fig. 1. — Evolution of the velocity and angular radius of the shock wave of CI Cam 1998 
in the model of the nova shell interaction with CSM. Left (a): the velocity in the model A 
(solid line) and in the model B (dashed line). Right (b): the angular radius in the model A 
(solid line) and in the model B (dashed line). Interferometric data (crosses) are taken from 
(Mioduszewski & Rupen 2004). 
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Fig. 2. — X-ray flux in 3-20 keV band in the model A (a) and model B (b) compared with 
the observational data (Filippova et al. 2008) (crosses). Insets show electron temperature 
in the model (line) and observations (Filippova et al. 2008). 
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Fig. 3. — Peak synchrotron luminosity of the expanding spherical shell vs. time of the peak 
at the frequency 2 GHz. Numbers at the lines indicate the velocity in units of 1000 km s _1 . 
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Fig. 4. — Flux evolution at 2 GHz for the expanding clumpy spherical shell. 1 - standard 
model (see text); 2 - density of the relativistic component is twice as smaller; 3 - filling 
factor is twice as smaller; 4 - cloud size is twice as smaller. 
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Fig. 5. — Flux evolution at frequences 2 GHz, 8 GHz, and 15 GHz in the model A (a) and 
model B (b) compared to data (Clark et al. 2000). 
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Fig. 6. — Evolution of X-ray and radio emission in the nonspherical model of CSM compared 
with data. Left (a): the same as in Fig. 2 but for nonspherical model; shown are the total 
flux (thick line) and contribution of shocks in polar cones (thin line). Right (b): the same 
as in Fig. 5 but for nonspherical model. The inset shows evolution of the angular radius (in 
mas) of the polar shock. 



